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Abstract 

A viscosity profile across the entire fluid outer core is found by interpolating between 
measured boundary values, using a diff"erential form of the Arrhenius law governing pressure 
and temperature dependence. The discovery that both the retrograde and prograde free core 
nutations are in free decay (Palmer and Smylie, 2005) allows direct measures of viscosity at 
the top of the outer core, while the reduction in the rotational splitting of the two equatorial 
translational modes of the inner core allows it to be measured at the bottom. We find 
2, 371 ± 1, 530 Pa • 5 at the top and 1.247 ± 0.035 x lO" Pa • 5 at the bottom. 

Following Brazhkin (1998) and Brazhkin and Lyapin (2000) who get 10^ Pa • s at the top, 
10^' Pa • 5 at the bottom, by an Arrhenius extrapolation of laboratory experiments, we use a 
differential form of the Arrhenius law to interpolate along the melting temperature curve to 
find a viscosity profile across the outer core. We find the variation to be closely log-linear 
between the measured boundary values. 

The close agreement of the boundary values of viscosity, found by Arrhenius extrapolation 
of laboratory experiments, with those found from the free core nutations, and the inner core 
translational modes, suggests that core flows are laminar and that the returned viscosities are 
measures of their molecular values. This would not be the case in the presence of the vigourous 
turbulent convection sometimes postulated by dynamo theorists. 

The local Ekman number is found to range from 10"^ at the bottom of the outer core to lO"^'' 
at the top. Except in the very lower part of the outer core, Ekman numbers are in the range 
lO^'* to 10~^, or below, in which the laminar flows of numerical dynamos and laboratory 
rotating fluids experiments occur. 



We find explicit expressions for the reciprocal Q's at both boundaries and for the viscous cou- 
pling torques between the outer core and shell, and between the outer and inner cores. For the 
high viscosity in the F-layer outside the ICB found from the reduction in the rotational splitting 
of the two equatorial translational modes of the inner core, the inner core is found to be tightly 
coupled to the outer core with a negligible contribution to dissipation in the free core nutation 
modes. 

"^F^iblished electronically in arXiv.org>physics >physics.geo-ph, Cornell University Library, Ithaca, N. Y., Febru- 
ary 2, 2008. 
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1 Introduction 



Properties of Earth's deep interior, such as its elasticity, density, pressure and gravity have tradition- 
ally been obtained through the inversion of seismic observations. While these have been important 
to our understanding of Earth's internal structure, the viscosity of the outer fluid core is crucial to 
our understanding of its dynamics and the generation of the geomagnetic field. Direct observations 
and limits on viscosity have traditionally been much larger than those values found on the basis 
of extrapolations of laboratory high pressure and temperature experiments (Lumb and Aldridge, 
1991). The latter tend to be close to that of liquid iron at atmospheric pressure (Rutter et al., 2002), 
while the former are many orders of magnitude larger (Davis and Whaler, 1997). 

Unusual properties of the lower outer core have long been suspected, going back to the 1926 
claim by Jeffreys (Jeffreys, 1926) of a strong negative P-wave velocity gradient there. While the P 
gradient is now thought to be small, but slightly positive, possibly due to solid inclusions slowing 
compressional waves in that region (Garland, 1971, pp.42-50), its properties remain a subject of 
speculation. A new seismic phase, PKhKP, has even been suggested by Bolt (BuUen and Bolt, 
1985, p.317) as originating from reflections in the lower outer core. 

With the discovery that both the retrograde and prograde free core nutations are in free decay 
(Palmer and Smylie, 2005), direct measures of the viscosity at the top of the outer core can be 
made. In this paper, we present a detailed analysis of both the Ekman layer at the top of the outer 
core and that at the bottom, just outside the inner core boundary. Our analysis yields a mean value 
of 2, 371 + 1, 530 Pa • s for the recovered dynamic viscosity at the top of the outer core. From the 
reduction in the rotational splitting of the two equatorial translational modes of the inner core, two 
independent measures of viscosity in the F-layer at the bottom of the outer core are found (Smylie, 
1999; Smylie and McMillan, 2000), as the reduction is larger for the retrograde mode than for the 
prograde mode. The retrograde equatorial mode gives 1.190±0.035xl0^^ Pa- 5, while the prograde 
equatorial mode gives 1.304 + 0.034x10^^ Pa- j for an average value of 1.247 + 0.035x10" Pa s. 

Our values of viscosity at the boundaries of the outer core are in close agreement with an 
Arrhenius extrapolation of laboratory experiments by Brazhkin (1998) and by Brazhkin and Lyapin 
(2000) who find 10^ Pa ■ s at the top of the outer core and 10" Pa - 5 at the bottom. Although 
the Arrhenius model is widely used over the limited ranges of pressure and temperature involved in 
laboratory experiments (Dobson, 2002), its exponential nature yields values of viscosity that appear 
extreme when extrapolated from laboratory pressures of the order of 30 GPa by more than a factor 
of ten to outer core pressures, and when the temperature extrapolation from laboratory values of the 
order of 2, 000 ^ is by a factor of two to outer core melting temperatures. Our measured viscosities 
at the boundaries appear to confirm the validity of Arrhenius extrapolation. 

The close agreement of the boundary values of viscosity with an Arrhenius extrapolation of 
laboratory experiments prompts us to construct a viscosity profile across the entire outer core using 
the Arrhenius law. Due to the strong pressure dependence of the activation volume, it is necessary 
to use a differential form of the Arrhenius law to interpolate along the melting temperature curve 
between the measured boundary values of viscosity. 

Both the free core nutations and the translational oscillations of the inner core are highly in- 
fluenced by Earth's rotation. They are examples of dynamical phenomena in contained rotating 
fluids. This is a subject that was of much interest several decades ago based on the inertial wave 
equation (Stewartson and Roberts, 1963; Roberts and Stewartson, 1965; Greenspan, 1969; Busse, 
1968; Aldridge and Toomre, 1969). In these descriptions, the container rotation is constrained 



2 



rather than free, there is no inner body, the fluid is assumed incompressible, uniform and non self- 
gravitating and the container is taken to be rigid. Although the theories were beautifully confirmed 
by laboratory experiments (Aldridge, 1967) their relevance to the real Earth in which the rotation of 
the container is unconstrained, there is an inner body, the fluid is compressible, stratified and self- 
gravitating and the boundaries are deformable does not appear to be direct. The governing equation 
then becomes the subseismic wave equation rather than the inertial wave equation (Smylie, Jiang, 
Brennan and Sato, 1992). The free core nutations are solutions of the subseismic wave equation 
which are close to pure rotations with respect to the shell, though not exactly so (Jiang, 1993; Jiang 
and Smylie, 1995; Jiang and Smylie, 1996). Thus, the motions in the Earth frame imitate the 'spin- 
over' mode of the inertial wave equation. Some of the basic features of the Ekman layers at the two 
boundaries in the real Earth resemble those of the 'spin-over' mode but previous theories are not 
given in detail with only results shown, possibly due to Greenspan's dictum that the 'computation 
of viscous eff'ects, is laborious but straightforward' (Greenspan, 1969, p.66). We give full details 
of the analysis of the Ekman layers at the two boundaries independent of the inertial wave equation, 
and calculate in detail the viscous coupling between the outer and inner cores and the shell in two 
Appendices. These details appear not to have been published previously. Explicit expressions are 
obtained for the reciprocal 2's at both boundaries, and for the viscous coupling torques between 
the outer core and the shell (mantle plus crust), and between the outer and inner cores. The inner 
core is found to be tightly coupled by viscosity to the outer core in the free core nutations. 

While the variational calculations of Jiang (1993) showed that there should be a prograde free 
core nutation (PFCN), in addition to the classical retrograde free core nutation (RFCN), he was 
anticipated in this discovery by both Mathews et al (1991), and by de Vries and Wahr (1991), in 
studies of the response of the Earth to nutational forcing. At the time, the VLBI nutation record 
was too short to give an acceptable confidence level, and no claim of observational support was 
made. 

There has been much interest in the detection of the inner core translational triplet of modes in 
the spectra of superconducting gravimeter observations since the initial identification fifteen years 
ago. They were first identified visually in the Product Spectrum of a total of 110,000 hours of 
observations at four European stations (see Figure 9 of Smylie, Jiang, Brennan and Sato, 1992). 
By adjustment of the degree one internal load Love numbers, they were quickly shown to have 
the correct splitting to be the translational triplet. A more complete study was made by Smylie, 
Hinderer, Richter and Ducarme (1993), including a complete analysis of the statistics of the Prod- 
uct Spectrum and a search across 4, 119 frequencies for correctly split triplets in the subtidal band 
between 2hr and Shr. The details of this search are summarized in Section 3 of the present paper. 
It was shown that the triplet originally identified visually had only one chance in 6.8 x 10^^ of being 
random. Nonetheless, the original identification generated a flurry of papers, pro and con, on this 
result. Hinderer, Crossley and Jensen (1995) were able to confirm the earlier result in the Product 
Spectrum of the same data sets, although by a different spectral analysis technique. Smylie, Hin- 
derer, Richter and Ducarme (1993) had used the Welch Overlapping Segment Analysis (WOSA) 
method, in which 12, 000 hr segments, windowed with a Parzen window with 75% overlap, en- 
abled an improved signal to noise ratio. In this case, the variance is inflated by a factor of 1.5 
but the record length is quadrupled, for an overall increase in eff'ective record length of 8/3. Hin- 
derer, Crossley and Jensen (1995) were unable to find the triplet in the cross spectrum of two years 
of common observations at just two stations, Cantley, Canada and Strasbourg, France. They had 
hoped for an improved result, even though simultaneous, short records from only two stations were 
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used, because the cross spectrum takes phase information into account. 

Apart from the difficulties associated with the low signal to noise ratio of the translational res- 
onances in the spectra of superconducting gravimeter observations, a second source of confusion 
is the calculation of the translational mode periods for a given Earth model. Most calculations use 
truncated vector spherical harmonic expansions derived from normal mode theory (Rogister, 2003). 
The problem with this approach is that the Coriolis acceleration couples these modes into slowly 
convergent infinite chains. Typically, these infinite chains are truncated after only a few terms. 
Johnson and Smylie (1977) showed that at periods of hours and more, convergence is very slow 
and that many terms need to be included for accurate computation of periods. This led Smylie, 
Jiang, Brennan and Sato (1992) to develop a variational, finite element method of computing long 
period core modes using bicubic splines as support functions. This method was used in the calcula- 
tion of translational mode periods and was extended by Jiang (1993) to determine solutions for the 
free core nutations. The finite element method, using local basis functions rather than global spher- 
ical vector harmonics, avoids the Coriolis coupling problem and allows the accurate calculation 
of periods. The initial formulation was criticized by Crossley, Rochester and Peng (1992) on the 
grounds that dynamical terms arising from pure translations had been neglected. This shortcoming 
was easily corrected (Smylie and Jiang, 1993). They showed that the degree one reciprocal Love 
number for the inner core has a quadratic dependence on frequency, as illustrated in Figure 5 of 
that paper. This shortcoming did not affect the original identification of the translational triplet, 
and current computations take the frequency dependence of Love numbers, used to describe the 
deformations of the shell (mantle and crust) and inner core, fully into account. The inviscid trans- 
lational mode periods for four Earth models are tabulated in Table 2, and are used in the splitting 
law plot shown in Figure 6. The computation of these periods is now done routinely to six or seven 
significant figures. 

Nonetheless, there remains a wide range of quoted translational triplet periods in the lit- 
erature. Crossley, Rochester and Peng (1992) give the triplet for Earth model 1066A as 
(4.04970, 4.43806, 4.89647) hours, close to the known values in Table 2 of (4.0491, 4.4199, 4.8603) 
hours. The next year, Crossley (1993) found (4.127,4.533,5.016) hours, while Rogister (2003), 
using normal mode theory, got (4.129,4.529,5.024) hours. Similarly, Rosat, Rogister, Crossley 
and Hinderer (2006) quote (4.64776, 5.24395, 5.83488) hours and (4.74744, 5.35865, 5.92142) for 
Earth model PREM, compared to the known triplet (4.6776, 5.1814, 5.7991) hours in Table 2. The 
noise level of quoted periods in the literature, even among papers with common authors, seems to 
exceed the noise level in the spectra of superconducting gravimeter observations by a large margin! 
While the scatter of reported periods may be due to attempts to apply truncated normal mode the- 
ory, the situation was complicated, even further, by Rieutord (2002) who finds (3.894, 4.255, 4.687) 
hours for Earth model 1066A. He misquotes the periods of Crossley, Rochester and Peng (1992) 
as (4.95, 4.438, 4.896) hours. The divergence of Rieutord's periods from known values, apart from 
typographical errors, appears to be due to his assumption that the shell (mantle and crust) and inner 
core are perfectly rigid, and that the fluid outer core is perfectly adiabatically stratified. He further 
dismisses the measurement of viscosity by Smylie and McMillan (2000), from the reduction of 
rotational splitting of the two equatorial translational modes, on the grounds that second order Ek- 
man boundary layer theory is required, even though he fails to perform the required second order 
analysis to support this claim. For the observed Ekman number of 1.2 x 10"^, since the Ekman 
layer theory is an ascending series in the square root of the Ekman number, second order theory 
would give a correction of 0.1 1 or 11%. The two independent measures of viscosity from the two 
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equatorial modes differ by 9.6%. It is not known if this difference is due to the neglect of second 
order terms in the Ekman boundary layer theory. 

The most comprehensive analysis of superconducting gravimeter observations appears to be 
that of Courtier et al (2000). A total of 294, 106 hours were assembled for analysis. In addition to 
Product Spectra, a global multistation experiment was performed using simultaneous observations 
from Brussels, Cantley, Kakioka, Strasbourg and Wuhan, allowing the use of phase to lower noise 
levels. Four independent sets of translational mode periods were recovered, including those shown 
in Figure 3 of the present paper, all differing from the original identification by no more than in the 
fifth significant figure. A similar global experiment using the Product Spectrum was attempted by 
Kroner, Jahr and Jentzsch (2004) for observations at five stations without result. Insufficient detail is 
given to determine if the data sets were correctly windowed to prevent frequency mixing from finite 
record effects, but the presence of tidal lines in the subtidal band suggests not. Also the inclusion 
of the Boulder data set, which has known serious timing errors, may have compromised the results. 
Finally, little is known about the excitation of the translational modes. Like the free core nutations, 
they may not be continuously excited, and may not be present in more recent observations. 

By contrast, a recent new development in the search for the translational modes has been 
achieved by Pagiatakis, Yin and El-Gelil (2007). Instead of using the Product Spectrum based 
on 12, 000 hour segments, across many stations, they use only the Cantley record of seventy-two 
months length. For each month they compute a least squares periodogram. Then, in a variation 
of the Product Spectrum method, they use the geometric mean of the seventy-two periodograms to 
look for translational modes. While the use of only one month data segments limits resolution, they 
find support for the periods found in the original detection using the Product Spectrum. 

As illustrated in Figure 7, the translational mode periods are very sensitive to inner core den- 
sity. From Table 3, we see that the Cal8 Earth model of Bolt and Uhrhammer gives periods very 
close to those observed. A reduction of only 2.25 milligrams • cm~^ in inner core density brings 
the three periods for this Earth model into coincidence with the observed triplet (Smylie, Francis 
and Merriam, 2001). In terms of the density jump at the inner core boundary, the Cal8 value of 
1, 170 kg • is reduced to 1, 168 kg • m~^. 

Traditionally, the dichotomy between extrapolated laboratory values of viscosity and directly 

observed values for the Earth's core (Lumb and Aldridge, 1991), the latter being many orders of 
magnitude larger, has been explained by dynamo theorists as due to turbulent flow, with the ob- 
served values reflecting eddy viscosities. It is difficult to see how turbulent flow could explain 
observed viscosities spanning nearly eight orders of magnitude across the core. This view is also 
challenged by the close agreement of our measured boundary values with the Arrhenius extrapola- 
tion of laboratory experiments by Brazhkin (1998) and by Brazhkin and Lyapin (2000). The latter 
reflect molecular viscosities. While turbulent flow may enhance dynamo action, successful numer- 
ical dynamos with laminar flow, and with Ekman numbers in the range implied by our measured 
viscosities have been routinely achieved (Olsen and Christiansen, 2002). 

2 Viscosity at the Top of the Outer Core 

In realistic Earth models, the finite element-based variational calculations of Jiang (1993) showed 
that there are two free core nutations. In addition to the classical retrograde mode (RFCN), a second 
prograde mode (PFCN) appears. The motions are illustrated by Poinsot constructions, in which a 
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small body cone rolls once per sidereal day without slipping on a large space cone as shown Figure 
1 below. 



Figure 1: Poinsot constructions for the Free Core Nutations. The RFCN is shown on the left, the 
PFCN on the right. In each mode, the small body cone rolls once per sidereal day without slipping 
on the large space cone. The line of contact is Earth's instantaneous rotation axis. 



In the Earth frame, both modes appear as nearly diurnal retrograde wobbles, the wobble asso- 
ciated with the RFCN is a little faster than retrograde diurnal, while the wobble associated with 
the PFCN is a little slower For the RFCN, the ratio of the nutation amplitude to the wobble 
amplitude a-^ is given by a^/aw - 1 = -^/ctn, where Q is Earth's angular rotation velocity and 
cr^v is the nutation angular frequency. For the PFCN, the ratio is given by a^/aw + 1 = CI/(Tn. 

In the study of Palmer and Smylie (Palmer and SmyUe, 2005), both the RFCN and the PFCN 
were found to be in free decay, on the basis of VLBI nutation measurement series from GSFC and 
the USNO in excess of twenty-three years length. The GSFC series runs from August 3, 1979 to 
March 6, 2003 and is taken to span a period of 8,617 days, while the USNO series runs from August 
3, 1979 to March 29, 2003 and is taken to span 8,631 days. In order to investigate the decays in 
detail, each record was divided into 2,000 day segments, advancing down the time axis in 400 day 
steps. Spectral densities were estimated on the basis of each of four successive 2,000 day segments. 
The result was spectral estimates centered at 1,600 days into the record, and at 400 day increments 
down the time axis thereafter The spectral amplitudes are shown plotted against time in Figure 2. 

In free decay, the logarithm of the nutation amplitude decreases linearly with time and 
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Figure 2: Amplitudes of the free core nutations plotted on a logarithmic scale as functions of time. 
Upper plots show the RFCN, upper left for GSFC, upper right for USNO, lower plots show the 
PFCN, lower left for GSFC, lower right for USNO. Linear fits to the time dependencies are shown 
directly on the plots, along with extrapolated values at the time origin. 



where is the signed nutation period (negative for retrograde, positive for prograde), Qn is the 
apparent Q of the nutation in the space frame, is the amplitude at time ? = 0, r is the e-folding 
time and ti/2 is the half life of the decay. The upper positive sign applies to the RFCN, while the 
lower negative sign applies to the PFCN. 

The actual physical dissipation takes place in the Earth frame, through the associated nearly 
diurnal retrograde wobbles, and is measured by Qw- Qw is related to the apparent nutation Q by 

Qw = ±(^-Tn/T,)Qm, (4) 
with Ts being the length of the sidereal day. The wobble Q is found to be 

Gw = -^^(l-r,™. (5) 

Results of the regression on the logarithm of the nutation amplitude, ([I]), are given in Table 1. 

In realistic Earth models, the two free core nutations have associated nearly diurnal retrograde 
wobbles that are close to pure rotations with respect to the shell (mantle and crust) and the inner 
core (Jiang, 1993; Jiang and Smylie, 1995; Jiang and Smylie, 1996). Their angular frequencies 
are close to -Q., the negative of Earth's angular rotation velocity. In spherical polar co-ordinates 
(r, 9, <p) the velocity field is closely 

V = -OAr sin (0 + Q.t) - ^Ar cos d cos (0 + Q.t) , (6) 
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Table 1 : Fitted decay parameters of the Free Core Nutations. 

c{lO^'^days~^^ d aNpiMas) TN{days) Qn Qw ti/2(days) 

RFCN 



PFCN 



GSFC 
USNO 

GSFC 
USNO 



-1.23909 
+0.1006 
-1.20360 
+0.1486 

-1.03217 
+0.1461 
-0.85237 
+0.1236 



2.69740 
2.74547 

2.51192 
2.37147 



498.20 
556.51 

325.03 
235.22 



-440.865 
+31.519 
-410.147 
+17.706 

474.308 
+86.578 
444.847 
+74.677 



26.404 
+4.031 
27.444 
+4.573 

34.052 
+11.036 

41.356 
+12.939 



11,068 
+2,481 
11,394 
+2, 390 

13,221 
+6,698 
16,009 
+7, 696 



2,429 
2,501 

2,916 
3,532 



for nearly diurnal wobble amplitude A. 

On the assumption that the boundary layers are a small fraction of the radius in thickness, 
the leading order boundary layer equations in the {6, (p) components of the extra boundary layer 
velocity, {ve,v^, are (Moore, 1978; Smylie and McMillan, 1998) 

dve d^Vg 
— -2Qv,cose=v— , 

-^ + 2Qv,cose=y-^, (7) 
at or^ 

with V denoting the kinematic viscosity. 

The detailed solution of these boundary layer equations, and the calculation of the rates of 
energy dissipation in each of the respective layers at the boundaries of the outer core is left to 
Appendix A. From equation (l85l) of that Appendix, the rate of energy dissipation in each of the 
respective boundary layers is found to be 



with po denoting the density, A the amplitude, and ro the radius at the boundaries. If is the 
amplitude of the nearly diurnal retrograde wobble of the outer core with respect to the inner core, 
and if Ah is the amplitude with respect to the shell, the total energy dissipated per cycle in both 
boundary layers is 

£ = • ^ = 33^' (po (a) Ay + Po ib)Ay V^) (9 V3 + 19) , (9) 

with Po (a) and representing the density and kinematic viscosity just outside the ICB (at ro = a) 
and Po (b) and Vh representing the density and kinematic viscosity just inside the CMB (at ro = b). 
The quality factor of the nearly diurnal wobbles accompanying the free core nutations, Qw, is 
defined as 2n times the ratio of the total energy to the energy dissipated per cycle, E. The total 
energy of the motion is closely 

^IocA\ (10) 
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where loc = 91 1.79 x 10^'* kg ■ is the moment of inertia of the outer core. The reciprocal of the 
overall quality factor is then the sum of the reciprocals of the effective quality factors at the two 
boundaries, 

^ = ^ + ^. (11) 

Qw Qa Qb 



with 

1 Inpo (a) V2i^ (9 VS + 1 9) 

Q'a^ 357;^^ ' 

1 27rpo(Z7)Z7^V2n7n(9V3 + 19) 



(12) 
(13) 



Neglecting the perturbation on shell rotation, if the amplitude of the nearly diurnal wobble of 
the outer core is B and that of the inner core is C, 

Ab = B, Aa = B-C. 

The detailed calculation of the viscous coupling torques the outer core exerts on the shell and 
inner core is left to Appendix B. From expressions (11031) and (11041 ) of that Appendix, and ignoring 
the small Chandler resonance effect, the equation of motion of the outer core is 

7„e-'"' (C - 5) - yte-'^^B = [b - iQB) e-'"", (14) 

while that of the inner core is 

- r,e-'"' (C-B) = /,-, (C - iQC) e-'"", (15) 

with lie representing the moment of inertia of the inner core. Equations (fT4l) and (fTSl) constitute a 
linear, homogeneous differential system. For time dependence e'^', it becomes 



A-ia + {ya + 7b)/Ioc -7a/ loc \( B 
-yjit, A-ia + rJIiJ\ C 



= 0, (16) 



with characteristic equation 



^ +1 — J ^Tr" 7~r" ' ^ ^ 



where A' = A - iQ,. 

From (fT03]) . (fT04l) . the ratio 

2^ _ po (b) M ^/n 
7a Po (a) Vv^ 



= 1/117.4, (18) 



for density ratio 0.8069, radius ratio 2.8668, and the ratio of the square roots of the kinematic 
viscosities 1/6400. Correct to terms of first order in the small quantity 7h/7a, the roots of the 
characteristic equation are 

A[ = ^, and A', = -^7a - ■^7b, (19) 
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with Ic = loc + he representing the moment of inertia of the entire core. The admissible solutions of 
the system (fT6l) are then the linear combinations 



Ce 



-iSlt 



J a Ic 



7h h 



J a Ic 



(20) 



with or, representing arbitrary linear combination coefficients. The decay times of the two solu- 
tions are, respectively 



Tl = 



Rljb 



loclic 



(IcRlya + licRlybY 



(21) 
(22) 



where 



V2 

Rlya =npo (a) a' ^^^^ [9 Vs + 19] , 
Rljh =npo (b) b' ^J^— [9 Vs + 19] . 



35 

Again, correct to first order in the small quantity jblja, the ratio of the decay times is 

T 2 lochc Rljb lochc lb 



Tl /2 Rlja II Ja 

With loc = 911.79 X 10^^ kg ■ m^, = 6.16 x lO^^ j^g . ^nd, hence /, = 917.95 x lO^^ kg 
we have ti/ti = 1/17,613. Thus, the second solution damps out rapidly, and we are left with 

7 a Ic 



(23) 
(24) 

(25) 



2 

• m , 



(26) 



and 



B-C = -^^-^B. 

Ja Ic 



(27) 



From expression (1851) . the total rate of energy dissipation at both boundaries then becomes 

v2n 



The energy dissipated per cycle is 



— = np^ (b)Alb' (9 V3 + 19) 



35 



l + Z^^ 



7 a \I, 



(28) 



Q dt 35 
The total energy of the motion is 

1 2 1 



(29) 



1 + libi 

7 a Ic 



n2 



1 

= -Ic 

2 " 



7 a \Ic 



(30) 
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to first order in the small quantity yblja- To the same order, the wobble quality factor is 

35/Ji + r,/r«(4//,)'l 

Qw = 2nT/E = ^ (31) 

2npo (b) M VSv^yn (9 V3 + 19) 

Finally, correct to first order, the viscosity recovered from the observed Qw of the nearly diurnal 
retrograde wobbles is 

n25f,a\l+2yt/7a(Iic/Icf] 

Vb = -. (32) 

Sn^plib) (9 yf3 + 19)' Q\ 

We see from this expression that the inner core is tightly coupled to the outer core motion, and that 
the correction for dissipation in the lower boundary layer, represented by the quantity in square 
brackets on the numerator, differs from unity by only 3.8 x 10"^ and may be neglected. 

Using expression (l32l) . and the values of listed in Table 1, we recover a viscosity of 3, 038 + 
1, 362 Pa • 5 for the RFCN from the GSFC series, 2, 866 + 1, 203 Pa • 5 from the USNO series, 
2, 129 + 2, 157 Pa - s for the PFCN from the GSFC series and 1,452 + 1, 396 Pa • 5 from the 
USNO series. The mean value of the recovered dynamic viscosity is then 2,371 + 1,530 Pa ■ s. 
For a density of 9.82 x 10^ kg ■ mT^ at the top of the core, the corresponding kinematic viscosity is 
0.2414 + 0.1558 m^^-i. 



3 Viscosity at the Bottom of the Outer Core 

Near the bottom of the outer core, the viscosity in the F-layer (Smylie, 1999; Smylie and McMillan, 
2000) can be found from the reduction in the rotational splitting of the two equatorial translational 
modes of the Inner Core. The translational modes are observed in the Product Spectrum of global 
networks of superconducting gravimeters (Smylie, Hinderer, Richter and Ducarme, 1993; Courtier 
et al., 2000). In Figure 3, we show the three translational mode resonances found in the Product 
Spectrum, based on observations at Bad Homburg (24,272 hours), Brussels (83,892 hours), Cantley 
(32,992 hours) and Strasbourg (78,504 hours). Both the prograde and axial mode resonances are 
well above the 95% C.I. and the retrograde mode is just below this level of significance (Fig. 8 (b). 
Courtier etal, 2000). 



A much more stringent test of significance arises from consideration of the pressure and viscous 
drags on the Inner Core (Smylie and McMillan, 2000) which leads to a splitting law of the form 



n,2,^^(ii-i=o. 



(33) 



where T is the period, Tq is the unsplit period, is the length of the sidereal day, and g" is a 
dimensionless viscous splitting parameter. For the axial mode, the viscous splitting parameter is 
related to the inviscid splitting parameter g' by 



\Mi-M' ^ 



(34) 
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Figure 3: Product Spectra of (from left to right) the retrograde equatorial, axial and prograde equa- 
torial translational modes of the Inner Core. The prograde equatorial mode is near the large solar 
heating tide feature 5 6 at exactly six cycles per solar day. The recovered central periods are, re- 
spectively 3.5822 + 0.0012 hr, 3.7656 + 0.0015 hr and 4.0150 + 0.0010 hr. 



and for the equatorial modes by 



1 - - 



+ 



Ekficr) 



S\Mj+/3 M[ + a 
a, /3 are coefficients of the pressure drag on the Inner Core given by 



(35) 



and 



a = M'j 



1 3 
_ + _. 



Ml + {a/by Ms 



2 '^Mo + Ms(l-{albf) 

13 Mi + (a/bf Ms 
4 ^Mo + Ms(l-{a/bf)l 



(36) 



(37) 



M[ is the mass of the Inner Core, Mq is the mass of the outer core, Ms is the mass of the shell, 
and Mj = 4/37ra^po(a) is the displaced mass, po(«) the density at the bottom of the outer core, 
cr = a)/2Q. is the dimensionless Coriolis frequency corresponding to angular frequency co. f" (cr), 
p (cr) are dimensionless functions of cr given by 



and 



r icT) = js [(a + If' + (cr- If'] - y [(cr + 1)^^^ - (cr - if']^ , 
r (cr) = 1+24 (±cr + D^^^ - 16 (±cr + 1)^^^ - ^ [i±cr - if - (±cr + if']^ , 



(38) 



(39) 



with the upper sign referring to the retrograde mode, for which cr is positive and the lower sign 
referring to the prograde mode, for which cr is negative. 

For three candidate periods, (retrograde), Tc (axial) and Tp (prograde), the splitting equation 
(13^ provides the corresponding values of the dimensionless viscous splitting parameter, 
and ^p, for a given value of Tq. Thus, the whole frequency axis can be searched for correctly split 
resonances. For a resonance centred on frequency /y, its form at neighbouring frequencies ft is 



1 + 42 [(ft - fj) /fj]' 



(40) 



12 



For record segments of 12,000 hour length. Product Spectral estimates Si are spaced at intervals 
of 1/12,000 cycles/hour along the frequency axis. In the subtidal band, between 2 hr and 8 hr 
period, there are 4, 501 spectral estimates. For twenty-five spectral estimates centred on frequency 
fj with Q = 100, the misfit of (|40l ) to spectral estimate Si is 

eij = AjVij -Si. (41) 

The error energy of the fit is 

7+12 

^ = Z 4- (42) 



Minimizing the error energy of the fit gives 



7+12 ;+12 



i=j-i2 i=j-n 



with minimum error energy 



7+12 j+n 



;=7-12 i=j-U 



As a measure of the strength of a potential resonance of the form (l40l).we use the parameter Sj = 
Ayi^in- When a large well-fit resonance is found, we expect 5^ to be large, and if a small, poorly 
fit spectral feature is found, we expect 5^ to be small. For each of the available 4,477 frequencies, 
fj, in the subtidal band, we set Tq = l/fj and compute = 1/Tr, fc = l/Tc and fp = -l/Tp 
from equation (|33l) . The values of 5^, 5^ and 5p of the resonance parameter 5^ at the discrete 
frequencies nearest fp, fc and fp, respectively, are then multiplied together to form the splitting 
product Pj as an indicator of the presence of correctly split resonances. We show the resulting 
probability density function (PDF) for the splitting product computed at 4, 119 points along the 
frequency axis in the subtidal band in Figure 4. 



500 




Splitting Product P 



Figure 4: Probability density function (PDF) for the splitting product P. Bins in P are 0.1 wide. 
The fitted PDF is for axl distribution with v = 2.97614 for the random variable 4.5314/'. 
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The PDF shown in Figure 4 allows the evaluation of the significance of translational triplets 
along the frequency axis. In Figure 5 we show the splitting products found between 2 hr and 10 hr 
periods. A very large value of P is found at Tq = 3.7975 hr. From the PDF it is found that the 
probability of a realization of Pj larger than the largest shown in Figure 5 is only 1 in 6.8 x 10^^! 
The resonances shown in Figure 3 seem to have been correctly identified as translational modes. 
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Figure 5: Splitting product P as a function of frequency. The large spike at Tq = 3.7975 corresponds 
to the translational triplet plotted in Figure3. 



n I I I I I r 



J I I I I I L 



Plots of the splitting law for the three translational modes are shown in Figure 6. The inviscid 
periods for the four Earth models plotted in this Figure are listed in Table 2. 

Table 2: Inviscid periods for the four Earth models shown plotted in Figure 6. 



Periods 


Retrograde 


Axial 


Prograde 




{hours) 


{hours) 


{hours) 


Corel 1 Periods 


5.1280 


5.7412 


6.5114 


PREM Periods 


4.6776 


5.1814 


5.7991 


1066 A Periods 


4.0491 


4.4199 


4.8603 


Cal8 Periods 


3.5168 


3.7926 


4.1118 



Two independent measures of viscosity are given, as the reduction in rotational splitting is larger 
for the retrograde mode than for the prograde equatorial mode. The retrograde equatorial mode 
gives 1.190+0.035x10" Pa -5, while the prograde equatorial mode gives 1.304+0.034x10" Pa-s. 
A balanced error value of 1.247 x 10" Pa • 5 yields viscous periods that are only 6.5 s longer than 
the observed periods. 

From Figure 6, we see that the observed periods are close to those for the Cal8 Earth model. In 
Table 3 a detailed comparison of the Cal8 periods and those observed is shown. 
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3.0 4.0 5.0 6.0 7.0 

Unsplit period Tq in hours 

Figure 6: Splitting curves for the three translational modes. The inviscid curves are shown dashed 
using the splitting parameters for Earth model Cal8 (open squares) of Bolt and Uhrhammer (BuUen 
and Bolt, 1985, Appendix). Inviscid periods are over plotted for Earth models Corel 1 (open stars) 
(Widmer et al., 1988), PREM (open diamonds) (Dziewonski and Anderson, 1981) and 1066A (open 
triangles) (Gilbert and Dziewonski, 1975). Solid viscous splitting curves are for a single viscosity 
of 1.247 X 10" Pa s. 

Table 3: Comparison of the observed translational mode periods with those of the Cal8 Earth model. 



Periods 


Retrograde 


Axial 


Prograde 




(hours) 


(hours) 


(hours) 


Observed Periods 


3.5822 


3.7656 


4.0150 


Cal8 Viscous Periods 


3.5840 


3.7731 


4.0168 


Cal8 Inviscid Periods 


3.5168 


3.7926 


4.1118 



The close match of the observed periods to those of the Cal8 Earth model is due to the sensitivity 
of the translational mode periods to inner core density. In Figure 7, we show the density profiles 
of the Inner Core for Earth models Cal8, 1066 A, PREM and Core 11 together with their unsplit 
periods, Tq. 



The axial mode period suffers little rotational or viscous splitting. Its observed period provides a 
strict constraint on inner core density. The calculated axial mode period for Cal8 is only 27 s longer 
than the observed period. An overall density decrease in the inner core of only 2.25 milligrams ■ 
cm~^ would bring them into coincidence, giving a very strong confirmation of Cal8. In Table 4 we 
show a comparison of the unsplit period, Tq, for the four Earth models with the observed value. 
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Figure 7: Detailed density profiles of the Inner Core for Earth models Cal8, 1066A, PREM and 
Corel 1. The range of 0.6 gm ■ cm~^ causes a nearly 2 hr difference in the unsplit period, Tq, giving 
a resolution of 200 minutes /gm • cm~^. 



Table 4: Comparison of the unsplit period, To, for four Earth models with the observed value. 



Earth Model 


Unsplit Period 


Deviation 


Error 




To (hrs.) 


ATo (hrs.) 


% 


Observed 


3.7985 






Cal8 


3.82470 


0.0262 


0.69 


1066 A 


4.45471 


0.65621 


17.28 


PREM 


5.23836 


1.43986 


37.91 


Corel 1 


5.81973 


2.02123 


53.21 



4 A Viscosity Profile for the Outer Core 

The boundary values of viscosity we have found are in very close agreement with an Arrhenius 
extrapolation of their laboratory experiments by Brazhkin (Brazhkin, 1998) and by Brazhkin and 
Lyapin (Brazhkin and Lyapin, 2000), who find 10^^ Pa ■ s at the bottom of the outer core and 
10^ Pa • 5 at the top. We are prompted by the very close agreement of the viscosity measures at the 
boundaries with those provided by the Arrhenius extrapolation, to extend it to interpolate between 
the boundary values, to obtain a viscosity profile across the entire outer liquid core. 

The Arrhenius description of the temperature and pressure dependence of the dynamic viscosity 
77 is (Brazhkin, 1998) 

)' ^^^^ 

with Eacto representing the activation energy at normal pressure, P the pressure, Vact the activation 
volume, k Boltzmann's constant and T the Kelvin temperature. Vact is proportional to the atomic 
volume, which, in turn, is inversely proportional to the density p. While the activation volume 
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for liquid metals at atmospheric pressure is very small, Brazhkin (1998) and Brazhkin and Lyapin 
(2000) report experimental results on pure iron at the melting temperature, r„,, that show it to be 
strongly pressure dependent up to pressures of 95 kbar. The strong pressure dependence requires 
integration of the differential form of the Arrhenius expression. For dominant pressure dependence, 
from expression (|45l) the differential increment in viscosity is proportional to 

D P 

expD—^dP, (46) 



pTm PT„ 

with D a pressure-dependent parameter, allowing for the pressure dependence of the activation 
volume, and dP being the differential increment in pressure. The integral of (|46|) over pressure 
is easily converted to an integral over radius r since dP/dr = -pg, where g is the gravitational 
acceleration at radius r. The viscosity at radius r is then 

r D IP \dP 
T](r) = T]h + T]h \ — expm— — Jr, (47) 
Jb pTm \ pTni/dr 

with b the radius of the core-mantle boundary and rjb = 2,2>1\ Pa ■ s the dynamic viscosity at the 
top of the core. To perform the integration in (1471) . we require profiles of pressure, density, melting 
temperature and pressure gradient. The pressure profile can be found by integrating the product 
of gravity and density for an Earth model (we use Cal8, see BuUen and Bolt (1985), p.472). The 
melting temperatures are found by spline interpolation onto the Cal8 radii from those tabulated by 
Stacey (1992, p.459). The required profiles are shown in Table 5. 



Table 5: Pressure, density, melting temperature and radial pressure gradient profiles. 



radius 


P 


P 


T 


dP/dr 


{km) 


(10" Pa) 


{10^ kg- m-3) 


(K) 


(lO^ Pa/m) 


1,216 


3.300 


12.20 


4,961 


-5.600 


1,371 


3.223 


12.14 


4,905 


-6.094 


1,571 


3.094 


12.03 


4,824 


-6.737 


1,821 


2.916 


11.84 


4,710 


-7.507 


2,171 


2.636 


11.52 


4,521 


-8.479 


2,571 


2.278 


11.11 


4,258 


-9.421 


2,971 


1.886 


10.62 


3,936 


-10.12 


3,171 


1.681 


10.33 


3,751 


-10.34 


3,371 


1.473 


10.01 


3,551 


-10.47 


3,486 


1.350 


9.860 


3,429 


-10.56 



Since the activation volume increases strongly with pressure (Brazhkin, 1998) as represented 
by our parameter D in equation (|46l) . we allow for a linear variation with depth through 

D = Dh + ^D,, (48) 
b - a 

where a is the radius of the inner core and Dh, are constants. The integration in (BTl ) is carried 
out by Simpson's rule over 100 steps with spline interpolation across the whole outer core. It is 
found that the constant, Db, controls the curvature of the viscosity profile near the core-mantle 
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boundary, while the curvature otherwise departs only slightly from log-linear. Some numerical 
experimentation shows that the profile is closely log-linear, even near the core-mantle boundary, 
for a value = 4.5 x lO""* m ^ ■ ■ K and that for = 2.976 x 10"^ m"^ ■ ■ K, the viscosity at 
the bottom of the outer core, rja = 1.247 X 10^^ Pa ■ s, is closely matched. The resulting viscosity 
profile is shown in Figure 8. 




n 



1216 1670 2124 2578 3032 3486 
Radius (km) 

Figure 8: Viscosity profile for Earth's Outer Core. 
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5 Discussion 



Given past uncertainties, the agreement between viscosities in the outer core measured from the 
VLBI observations of nutations, the superconducting gravimeter observations of the translational 
modes, and the Arrhenius extrapolation of laboratory high pressure and temperature experiments, is 
quite remarkable. The very large gap between direct observations and extrapolations of laboratory 
values appears to have closed. The viscosity values involved are for molecular viscosities and they 
are large enough that flows in the core are likely laminar in contradiction to conventional thought 
that the flows are turbulent and that the large viscosities reflect eddy viscosities. 

According to the viscosity profile we have derived, the local Ekman number ranges from 1.2 x 
10"^ at the bottom to 2.7 x 10" at the top of the outer core. Our results appear to confirm the 
suggestion by Braginsky (Braginsky, 1963) that the release of the latent heat of fusion as metallic 
constituents freeze out in the F-layer at the bottom of the outer core may be the energy source 
required to drive the geodynamo through compositional convection as studied by Loper and Roberts 
(Loper and Roberts, 198 1). At the same time, except in the very lower part of the outer core, Ekman 
numbers are in the range 10""* to 10"^ ,or below, in which numerical dynamos operate (Olsen and 
Christiansen, 2002). 
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Appendix A 



Ekman Boundary Layers and Dissipation 

In this Appendix, we describe in detail the solution of the boundary layer equations ([7]) and 
the calculation of the rates of energy dissipation in the two boundary layers at the boundaries of 
the outer core. For convenience in keeping track of phase, the velocity components (ve, v^pj will be 
taken to be complex phasors with time and longitude variations given by e~'^'^'^^'\ For an assumed 
radial dependence proportional to e'^'', we are lead to the homogeneous system of equations 



-2Qcos0 A^v + iQ. 



V0 



= 0. (49) 



For this system to have a solution, A must satisfy 

[a\ + inf = -4Q2 cos^ e. (50) 

Thus, 

A\ + in = +i2n cose. (51) 
Substitution of relation (|5T| ) into the system of equations (H^ . gives 

= +ive, (52) 

and 

= ±iv^. (53) 
The boundary layers are characterized by the dimensionless Ekman number 



(54) 



with length scale fixed by the radius b of the CMB. The four roots of the secular equation (|50l) are 
+Ai , +A2 with 



1-z /1/2 + COS0 1+?' ^ ,^ 

^i=^J^-^ = ^ for 0<2n/3 (55) 



1 + j /-1/2-COS0 1+?' ^ ,^ 

' ^-^ = — for e>2nl3 (56) 



and 



1-z /1/2-cose 1-z 

A2=—-J^— = — for 6>nl3 (57) 

b V Ek 62 



1 + z /-1/2 + COS0 1 + z 

' ^ = — for e<nll>, (58) 

Ek h 

with 61,62, the respective boundary layer thicknesses, both of O ( V^T)- 

In the boundary layer near the top of the core, the perturbing velocity components ^vg, vanish 
with decreasing radius and increase with radius, so that at r = Z?, they are equal and opposite to 
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the interior nutation velocity components ^ to satisfy the no-slip condition at the CMB. Thus, the 
admissible values of A satisfying equation (|5TI) are Ai,A2. In the boundary layer near the bottom 
of the outer core, the perturbing velocity components vanish with increasing radius and increase 
with decreasing radius to satisfy the no-slip condition at the ICB. There, the admissible values of 
A satisfying equation (|5T1) are -ii, -/I2. The perturbing velocity components are then given by the 
linear combinations 



(59) 
(60) 



where the lower signs apply for the range n/3 < 6 < InlZ, while the upper signs apply for the 
range < < 7r/3 for terms involving 82, and the upper signs apply for the range 27r/3 < 6 < n 
for terms involving 5\. Ar is the increment in radius. At the top of the core Ar = - r, and at the 
bottom of the outer core Ar = r - a. In general, the linear combination coefficients are complex 
with real and imaginary parts expressed by 



f=a + i/3 
g =7 + i^- 



(61) 
(62) 



Both the real and imaginary parts of expressions (|59l ) and (|60l ) are solutions of the boundary 
layer equations (|7]). Our interest is in the real parts of the velocity components given by 



-Ar/S, 



I Ar . Ar\ ^ 
a cos — + p sm — cos (0 -I- Qt) 
\ ^1 Si I 



, Ar . Ar . . ^ 

+ p cos h a sm — sm (0 -1- Qt) 

61 61 



+e 



-Ar/62 



Ar . Ar\ ^ ^ ^ 
y cos — + 6 sm — cos (0 -I- \lt) 

62 02 



. Ar Ar . 

+ I e cos h 7 sm — | sm (0 -I- Q.t) 

62 62 



(63) 



Rlv, 



-ArlSi 



Ar . Ar\ ^ 
-p cos — + a sm — | cos (0 -I- ilt) 

5\ 5\ 



. Ar . Ar . . ^ 
+ a cos — + p sm — sm (0 -1- ilt) 

61 61 



+e 



-ArlSi 



Ar . Ar\ ^ ^ ^ 
6 cos — -I- 7 sm — cos (0 -I- lit) 

02 02 



. Ar Ar\ 

+ I -7 cos 1- e sm — sm (0 -I- Q.t) 

62 62 



m 



For the velocity components (1631) and (|64|) to cancel the components of Q at the boundaries, we 
have 



Rlve ={a + y) cos (cp + Q.t) + (J3 + e) sin (0 + Q.t) = Aro sin (0 -I- Q.t) , (65) 
Rlv^ ={6-/3) cos (0 -I- Q.t) + ia-y) sin (0 + Q.t) = Aro cos 6 cos (0 -I- Qt) , (66) 
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where tq is the boundary radius, ro = b at the top of the core and ro = a at the bottom of the outer 
core. Then, 

/3 = ^Aro (1 - cos 61) , e= ^^^o (1 + cos 6) , a = y = 0. (67) 

Differentiation of expressions (|63l) and (|64l) yields the derivatives of the velocity components at the 
top of the core, 



dRlvg 

dr 
dRlv^ 

dr 



5i 52 



+— + — \cos{(f> + Qt) + \ — + — \sm((p + Qt) 



= -— + — cos (0 + Q.t) + +— + — sin (0 + Q.t) , 



6i 62 



(68) 
(69) 



while those at the bottom of the outer core are the negatives of these. At the top of the core, the 
leading order stresses on the outer surface are 



dRlva 



o-re = ri- 



dr 



and (T,-A = rj- 



dRlv. 



dr ' 



(70) 



where 77 is the dynamic viscosity. At the bottom of the outer core, the leading order stresses on the 
inner surface are the negatives of these, so that on both surfaces the rate of dissipation of energy 
per unit area in the motion with velocity I® against these stresses is 



de 

— =V0(Tre + V^CTr. 



=Aror] 
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+— + — cos ((/> + Q.t) sin (0 + + — + — sin (0 + 0.t) 



5\ 



+ cos^ll — - — |cos (0 + Q.t) + (+— + — |cos(0 + Q.t) sin(0 + Q.t) 
Substitution from ([67]), ([55]), ([561), dST]) and dSS]) yields, for < 6* < nl3. 



(71) 



^ =\poA^rl [ VvQ(l/2 + cose) 



{sin^ 6 (cos {(p + Q.t) sin (0 + Q.t) - cos^ (0 + Qf)) + 1 - cos 6\ 
-^jvQ.{-\|2 + cose) 

{sin^ e (cos {(p + Q.t) sin {(p + Q.t) + cos^ {(p + Qr)) - 1 - cos 0}] , 
and for 7r/3 <6 < In 13, 

^ =^PoAV2 [ VvQ(l/2 + cose) 

jsin^ (cos (0 + Qf) sin (0 + Q.t) - cos^ (0 + QO) + 1 - cos 
+ Vv^^(l/2-cos0) 

{sin^ e (cos (0 + sin {(p + Q?) - cos^ {(p + QO) + 1 + cos 0}] , 



(72) 



(73) 
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and for 2;r/3 < 6 < n. 



jsin^ (cos {(p + sin (0 + Q?) + cos^ {(p + fir)) - 1 + cos 
+ Vv^^(l/2-cose) 

[sin^ (cos + Q.t) sin (0 + Q.t) - cos^ + QO) + 1 + cos e^j] , 



with po representing the density just inside the respective boundary. 

Integrating over the entire boundary, the total rate of energy dissipation is 



dE 2 r^"" r de . 
—r =f'n — sin 9d6d(h 

dt "Jo Jo dt 



X 



+ 



+ 



+ 



+ 



J y^:^+<^o^s(^^'i'^0(\-co^e + %m^efjde 

I y-- + cos6'(sin0(l +COS0 + sin^e)) 
I J- +cos6'(sin6'(l -cose + sin^e)) ^6* 

Jn/3 ' ^ 

I J- -cos6'(sin6'(l +cose + sin^e)) J0 

^ ^-i -cos6'(sine(l -cos^ + sin^e)) 



J0 
^6/ 



+ J - cose(sin0(l + COS0 + sin^^jji 

where, in the integration over 0, we have made use of the integrals 

COS ((f) + Q.t) sin (0 + D.t)d(f> = - I sin {2 (0 + Q.t)} d(p = 0, 
2 Jo 



and 



cos^ (cf) + Qt) d(f> = - I (cos{2(0 + riO} + = TT. 
2 Jo 

We may write (1751) in the shorthand 



— = --poA'rl VvQ [7i + 72 + ^3 + 74 + ^5 + 76] , 



with 7i, 72, 73, 74, 75, 76 representing the six integrals in ( [75l ). Evaluation of the integrals 
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sion (1751 ) depends on the indefinite integrals 
/i = J" (I +cos0j ' sin^ede 

lis \3/2 , 8 (S \5/2 16 IS 

= ±-^-+cos0j sin^e- — ^- + cos0j COS0+— ^- + cos0j , (79) 

/2 = J" (I + cos ^ (1 - cos 0) sin OdO 
= ±-^-+cos0j +-^- + cos0j cos^-— ^-+cos0j , (80) 

h = ^ (^+^^^^]' (1 + '^os ^) sill 

2 /5 \3/2 2 /5 \3/2 4 IS \5/2 

"~3\2^'^°^^) ~3l2^^"^^) cose + + cos 0j , (81) 

s representing the sign, which can take on either the value + 1 , or the value - 1 , throughout each 
expression. Using (1791 ). (f80l) and (fSTT ). and inserting limits of integration, we find 

13 16 18 /3 57 r- 1 44 

•^^"30^T05"35V2' -^^-"105^^' " "6 " 105' ^^^^ 

1 44 , , 57 r- ^ ^ 13 16 18 /J , 

74 = = 73, 75 = V2 = 72 , 76 = — + J- = Ji. (83) 

^ 6 105 ^ ^ 105 ^ 30 105 35 V 2 ^ ^ ^ 

The six integrals sum to 

7i + 72 + 73 + 74 + 75 + 76 = -2 V2(9 V3 + 19) /35. (84) 

Substituting this result in expression (1781) . we find the rate of energy dissipation in each of the 
respective boundary layers to be 

^ = ^'^o^'^o Vm(9 V3 + 19) . (85) 
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Appendix B 



Viscous Coupling to the Inner Core and Shell 

From relation (fTTl) . it is apparent that the reciprocal of the overall quality factor cannot be 
less than the reciprocal of the effective quality factor arising from the boundary layer at the ICB. 
Because of the high viscosity there, the inner core is likely to be tightly coupled to outer core 
wobble. The nearly retrograde diurnal wobbles of the outer core, associated with the free core 
nutations, give rise to viscous torques exerted by the outer core on the inner core and shell. 

The viscous torques exerted by the outer core at the boundaries are 



J"r X (ecTre + ^(Tr,^ = [-6r(T,-^ + ^rcTre) dS, 



(86) 



where the integral is over the respective boundary surface. The spherical polar unit vectors are 

I, j, kj by 



6 = ICOS0COS0 + 7cos0sin0 - ksinO, 

4> = -i sin (p + j cos 0. 
The Cartesian components of the viscous torques are then 

-'In /-'TT 

r =ir, 



(87) 
(88) 



I sin ^f-cr^^ cos cos - cr^e sin 
Jo 

I sin Q \-(Tr^ cos sin + cr^e cos m d6d(p 
Jo 

Jr»2;r r^n 
I cr,-0 sin^ 6d6d(p. 
Jo 



+krl 



(89) 



From the expressions (|70l) for the leading order stresses, the integrations over (p are seen to depend 
on the elementary integrals 

Jr>2ff pin 
cos((p + Qt)cos^d^ = ncosQt, I sin + Qt) cos ^d(f) = n sin Q.t, 
Jo 

Jr>2ff p2n 
COS {(p + Q.t) sin (pd4> = -n sin Q.t, I sin {(p + Qt) sin (f)d(f) = n cos Q.t, 
Jo 

cos {(p + Q.f) d(p = I sin (cp + Q.t) d(p = 0. 
Jo 

The latter two integrals ensure that the viscous torques have only equatorial Cartesian components, 
(Yx, Yy), and writing f = F^c + iYy, we have 



Jo 



/3 



cos 6 sin 6 \ sin 1 



Jo 



COS sin +— + — + sin +— + — 



de 
de. 



(90) 
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Substitution from equations (|55] ). (I56l) . (1571 ). (|58l ) and (|67l) shows that evaluation of the torque 
expression (|90l ) depends on the integrals 



/V(F 



cos 6* I cos 6* sin 

\3/2 4 / ^ 
+ - COS 0( -T + COS 01 +— I-+COS0 



15 V2 



5/2 



cos^lcos^ sin 6d6 



3 \2 



15 \2 



\5/2 



2/^ .\^^^ 2. 8 ^, 

= +- -+COSe COS 0+—-+ COS COS0 + 



3/2 



16/5 \^/2 
(-TCOS^j 



105 \2 



+ cos^lsin^J^ ~ ~^ ^^^^ 



(91) 



where, again, 5 represents the sign, which can take on either the value +1, or the value -1, through- 
out each expression. 

The total viscous torque is made up of contributions from three zones of latitude. The contribu- 
tion from the region < 9 < n/3 is 



ro<e<;r/3 = - ^PoAro VvQe 



41 



-iSlt 



41 



3 19 1 



140 35 V 2 35 V2 



3 19 1 



140 35 V 2 35 V2 
while the contribution from the region n/3 < 6 < In 13 is 



Tni3<e<2ni3 = -np^Arl VvQe"'"'-- (-1 + /) , 



41 



70 



and that from the region 2nl3 < 6 < n h 



-iat 



41 



3 19 1 



140 35 V 2 35 V2 



+ — a/t 



3 19 1 



140 35 V 2 35 V2 



The total viscous torque is then 
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f = npoAri Vvae-'"'— V2 



35 



1 + / + — V3(l-0 



(92) 



(93) 



(94) 



(95) 



Separating the torques from the two boundary layers, the outer core exerts the viscous torque 



fa = npo {a)Aaa' V^e"'"'^^ <2 



\+i + —yf3{\-i) 



(96) 
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on the inner core, while it exerts the viscous torque 



19 



+ — V3(l-/) 



(97) 



on the shell. 

The extra nearly diurnal retrograde wobble of the outer core compared to its boundaries, in 
complex phasor notation, is 



oj = oj\ + ia>2 = Ae 



-iClt 



(98) 



with wobble angular velocity oj = (0)1,0)2). The rate at which the outer core does work against the 
viscous torques is 

YiOJi + T2OJ2 = ]r [ojf* + oft) , (99) 



dE 

— = T -OJ 

dt 



where the torque vector is F = (ri,r2), and where the superscript asterisk indicates the complex 
conjugate. From the torque expression (|95] ), we find the rate of energy dissipation in the boundary 
layers to be 

j7 = ^'^o^''"o V2^(9 V3 + 19) , (100) 

in agreement with equation (|85l) . 

The expressions (|96| ) and (|97] ) for the viscous torques the outer core exerts on its boundaries 
may be abbreviated to 

(101) 



fa = Yae'^^Aa 



and 
where 

and 



fb = J be '^*Ab, 
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ja = npQ (a) V^^^^^ 



■19 



1+/ + — V3(l-0 



1 + z + — V3(l -/) 
19 



(102) 
(103) 

(104) 
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